Group Velocity of Discrete-Time Quantum Walks 
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We show that certain types of quantum walks can be modeled as waves that propagate in a 
medium with phase and group velocities that are explicitly calculable. Since the group and phase 
velocities indicate how fast wave packets can propagate causally, we propose the use of these wave 
velocities in a new definition for the hitting time of quantum walks. The new definition of hitting 
time has the advantage that it requires neither the specification of a walker's initial condition nor 
of an arrival probability threshold. We give full details for the case of quantum walks on the Cayley 
graphs of Abelian groups. This includes the special cases of quantum walks on the line and on 
hypercubes. 
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I. INTRODUCTION 

While classical random walks have found important 
applications in classical computing, quantum random 
walks are anticipated to lead to significant applications in 
quantum computing. We will here focus on discrete-time 
quantum walks. These were first introduced by Aharonov 
et al [l[ as an example of a quantum physical process 
that exhibit striking differences from classical processes. 
Such differences have motivated a search for quantum al- 
gorithms based on quantum walks that outperform their 
classical counterparts @, H, 0| • 

Aharonov et al @ developed the basic theory of quan- 
tum walks on graphs with a focus on mixing-time prop- 
erties. They proved that quantum walks mix at most 
polynomially faster than classical random walks. Mar- 
quezino et al @ presented an analytical expression for 
the mixing time of discrete-time quantum walks on the 
hypercube. 

Particularly important in this context is the concept 
of hitting time. For classical random walks, the hitting 
time is unambiguously defined as the average time the 
walker takes to hit the final vertex for the first time af- 
ter departing from the initial vertex. The generalization 
for quantum walks is not straightforward, however, since 
measurements disturb the movement of the walker, see 
e.g. 0- One possibility is to let the walk evolve unmea- 
sured, i.e., unitarily, until the arrival probability at the 
final vertex is above some threshold. Another possibility 
is to perform a partial measurement at each step of the 
walk to check whether the walker has already reached 
the final vertex. Both definitions of arrival times have 
drawbacks. In the first case, the walker is not confirmed 
to have hit the final vertex, as there is only a probability 
of hitting it. In the second case, the quantum walk has 
been modified by the repeated measurements so that one 
is actually calculating the hitting time of a non-unitary 
walk that is effectively subject to a quantum Zeno effect. 
Nevertheless, at least for the symmetric walk on the hy- 



percube @, the two strategies for defining the hitting 
time yield similar results. We note that the hitting time 
can, in general, be infinite. Krovi and Brun [f| analyzed 
the conditions for infiniteness of the hitting time for walks 
on Cayley graphs based on the second definition. 

Our strategy here for defining the walker's speed is to 
view the quantum walker as a wave packet which evolves 
according to the Schrodinger wave equation and then to 
calculate the group velocity. In the case of the quantum 
walk on the line, for example, the problem is translation 
invariant. Therefore, the wave equation can be decom- 
posed into normal modes labeled by wave numbers, i.e., 
by wavelengths. More generally, one can consider, for ex- 
ample, quantum walks on graphs with a general Abelian 
symmetry group. There is then a corresponding gener- 
alized Fourier transform which leads to a normal mode 
decomposition labeled by generalized wave numbers. 

Using the normal mode decomposition, the unitary 
time evolution operator for the basic time step can then 
be diagonalized as a function of the wave numbers. 
The logarithm of the time evolution operator yields the 
Hamiltonian, and thus the energy as a function of the 
wave number. This is what we are after. The energy, 
divided by, or differentiated with respect to, the wave 
number yields the phase and group velocities, respec- 
tively. In this sense, quantum walks behave like waves 
propagating in a medium. For unitary quantum walks 
the waves are neither absorbed nor amplified, i.e., the 
propagation medium is passive. When the group veloc- 
ity is smaller than the phase velocity, the group velocity 
should then indicate the signal velocity, i.e., the velocity 
with which information can propagate. 

We can now define the notion of group velocity based 
hitting time for a quantum walk on a graph as the length 
of the minimal path connecting the initial and final ver- 
tices divided by the maximal group velocity for all wave 
numbers, assuming, as we will find is the case, that the 
group velocity remains below the phase velocity. We will 
compare the group velocity based hitting time with the 
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standard definitions of hitting times for quantum walks 
on one-dimensional lattices and on the rt-dimcnsional hy- 
percube. For the lattice, the hitting times coincide ap- 
proximately. For the hypercube, the group-velocity hit- 
ting time is 0{n^Jn), which is greater than the result 
from the standard definitions that is essentially 0(n) Q. 
We also analyze the group velocity of quantum walks on 
Cayley graphs of finite Abelian groups. 

The paper is organized as follows. In Sec. |TT] we re- 
view the theory of quantum walks on Cayley graphs of 
finite Abelian groups and present a general procedure 
for calculating the group velocity. In Sec. IIIII we present 
definitions of quantum hitting times. In Sec. HVI we com- 
pare the group-velocity hitting time with the standard 
definitions of hitting time for a quantum walk on the 
hypercube, and in Sec. [V] for a quantum walk on one- 
dimensional lattices. In the last Section we draw conclu- 
sions. 



II. QUANTUM WALKS ON CAYLEY GRAPHS 

A Cayley graph encodes the structure of a discrete 
group. Let G be a finite group and S a generating set. 
The Cayley graph T(G, S) is a directed graph such that 
the vertex set is identified with G and the edge set con- 
sists of pairs of the form (g, gh) for all g € G and h € S. 
The Cayley graph depends in an essential way on the 
generating set. It is interesting to diminish that depen- 
dence by demanding that S be a symmetric set, that is, 
if h G S then h^ 1 S S, where h~ Y is the inverse of h. 
In that case, the Cayley graph is an undirected regular 
graph of degree \S\ with no loops, where \S\ is the car- 
dinality of S. From now on we consider only symmetric 
generating sets. 

Coined quantum walks can be defined on T(G, S) in 
the following way. Let Hs be the Hilbert space spanned 
by states \h) where h G S. Hs is the coin or internal 
space. Let Hg be the Hilbert space spanned by states 
\g) where g G G. Hg is the physical stage where the walk 
takes place. The evolution operator for one step of the 
walk is U = S o (C ® I) where 



c= J2 c hihM{h 2 



(i) 



is the coin operator, / is the identity operator in Hg, and 
S is the shift operator given by 



S\h)\g) = \h)\gh). 



(2) 



As we see from the last equation, if the walker is in vertex 
g and the result of the coin toss is h, then the walker 
moves to its neighboring vertex gh. 

The analysis of the evolution is simplified in the Fourier 
space. Let us suppose that G is an Abelian group. In 
that case, G is a direct sum of cyclic groups, that is, 
there are integers Ni to N n such that G is isomorphic to 



Zn 1 x • ■ ■ x Zjv n , where is the additive group mod- 
ulo N. Any element g € G can be written as a n-tuple 
(ffir • " >fin)- Such decomposition can be determined ef- 
ficiently [9J. The Fourier transform on G is given by the 
operator 



F G 



I 
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where Xg is a character of G given by 



x 9 (h) = n 



(3) 



(4) 



where ui^. — exp(^P-) is the N 3 ; -primitive root of unity. 

The Fourier basis is an orthonormal set of vectors de- 
fined by 
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\Xh) = 
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(5) 



where h G G. That basis is interesting because any vector 
|/i)|x g ) is an eigenvector of the shift operator, in fact 



S\h)\x g )=Xh(g)\h)\x g ) 



(6) 



which can be proved by using Eq. ([J). If we analyze 
the form of the evolution operator, we conclude that in 
the Fourier basis it acts non-trivially only on the coin 
subspace. So, let us proceed with a reduced version of 
the evolution operator U g that acts on states \^f g (t)) = 
(x g \^{t)), where |^(£)} is the generic state of the walk 
at time t given by 



!*(*)) = EE th, g m)\x g ) 



(7) 



heS g£G 



U g acts on states \^ g (t)} = ^2 he s i'h^ityh) and the ma- 
trix components are {hi\U g \h 2 ) = Xhi(g)Chih 2 - 

The evolution equation is given by \^ g (t + 1)) = 
Ug\^ g (t)), which can be solved recursively yielding 
\V g (t)) = (Z7 fl )*|# fl (0)). One may calculate {U g f by di- 

agonalizing the matrix U g . Let us call exp(iwg' 1 ' 1 ), h G S 
the eigenvalues of U g . We can calculate the Hamiltonian 
H g using that U g — e~xjp(—iH g ) after taking h = 1. In 

the eigenbasis of U g , we have H g — diagjw^}. 

In Eq. (O, parameters h and g play a dual role. Pa- 
rameter g plays the role of the spatial position while pa- 
rameter h can be considered a generalized wave number. 
The differentiation of the energy lo^ with respect to h, 
or some quantity directly obtained from h, defines the 
group velocity. The ratio of the energy to h defines the 
phase velocity. 



III. HITTING TIMES 



Using the group velocity, v g , we can calculate the trav- 
eling time from vertex g\ to g 2 - Taking edges of length 
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one, the time is given by d/v g , where d is the length of 
the shortest path connecting the vertices. We define the 
group -velocity hitting time as the length of the shortest 
path divided by the maximal value of the group velocity. 

It is interesting to compare that "physical" hitting time 
notion with the mathematical definition that generalizes 
the well known classical hitting time notion. In the classi- 
cal case, the evolution is governed by a stochastic matrix 
and the hitting time is the expected time the walker takes 
to hit vertex g2 for the first time starting from vertex 171 . 
In the quantum case, there is more than one notion of 
quantum hitting time 0, El • Either one lets the walk 
evolve unitarily after leaving from vertex gi and checks 
when the probability at vertex gi is above some thresh- 
old, or one performs a partial measurement at each step 
to measure when the walker has reached vertex <? 2 - The 
first notion has the following definition. 

Definition III.l (One-shot hitting time) Given a 
threshold < p < 1 and a initial condition \(f>o) for the 
coin state, the one-shot hitting time from vertex gi to gi 
of the discrete-time quantum walk U is 

min{T| WeH s : ]T |(/i|(fl 2 |tf T |^o)|<7i)| 2 > P}- 
hes 

The hitting time may be infinite if one chooses p too 
high. On the other hand, it is advisable to take p as high 
as possible to have a good chance to find the walker on 
vertex gi. 

The second notion has two definitions. Let us first 
define the concurrent hitting time 0]. 

Definition III. 2 (Concurrent hitting time) A 

discrete-time quantum walk U has a (T,p) concurrent 
hitting time from vertex g\ to g2, if the \gi) -measured 
walk from U with the initial state \4>o)\gi) has a prob- 
ability greater or equal to p of stopping at a time 
t < T. 

A walk is called \g) -measured when we perform a mea- 
surement at each step of the evolution with the projec- 
tors P = I ® \g)(g\ and Q = I — P. If P is measured the 
process stops, otherwise the iteration is continued. 

Krovi and Brun [Io| proposed an alternative definition 
which does not have a threshold p. 

Definition III. 3 (Average hitting time) A discrete- 
time quantum walk U with initial state pi — \^){^\ where 
1^) = I0o)|<?i) has a (51,32) average hitting time 

00 

*=i 

where 

p(t) = T,{PU{QU) t - 1 Pl {U^Q) t - 1 U^P}, 
P = I®\92){g2\ andQ = I-P. 



Note that the wave function is not renormalized after the 
measurement at each step. 

A drawback of both the one-shot and the concurrent 
hitting times is that they depend on a choice of threshold 
probability. Intuitively the threshold should not be cho- 
sen too low, because else the hitting times would reflect 
the arrival of mere traces of probability. Exponentially 
suppressed traces of probability often arrive quickly but 
in practice cannot be considered a useful criterion for 
the arrival of the walker. The threshold probability also 
should not be set too high, as this could lead to an in- 
finite hitting time. Apart from these arguments it ap- 
pears difficult, however, to further constrain any choice 
of threshold probability. 

Let us consider, therefore, that, intuitively, the walker 
carries information and that it is necessary to wait for 
this information to arrive at the final vertex. Informa- 
tion travels in a medium with what is called the signal 
velocity which, for normal dispersion, should be given by 
the maximum value of the group velocity. This is the 
case here, as the analysis for walks on the line and on the 
hypercube will show. The medium is passive, since the 
evolution is unitary, and the dispersion relations are well 
behaved with the group velocity staying below the phase 
velocity. 

IV. QUANTUM WALKS ON HYPERCUBES 

The n-dimensional hypercube is the Cayley graph of 
the group Z^. Let us represent the group elements by 
binary n-tuples x — (x„_i, . . . x\, xq) and the generating 
set by {ej, < j < n}, where ej has a single 1 entry 
in the (n — j) th component. In this case the vertices 
(0, • • • ,0) and (1, • • • ,1) are opposite corners. The shift 
operator Q reduces to 

S\e j )\x) = \e j )\x@e j ), (8) 

where © is the n-tuple binary sum. The character is 
given by Xxief) = (— l) Xj and the matrix elements of the 
reduced evolution operator are (ei\Uk\ej) = {—l) ki Cij. 

From this point on, let us particularize the analysis to 
the n-dimensional Grover coin, C'ij = 2/n — Sij, which 
obeys the permutation symmetry of the hypercube. For 
this coin, we can calculate explicitly the eigenvalues of 
C/fe. They are given in the following table [6(. 



Hamming weight 


Eigenvalue 


1*1=0 


-1 
1 




-1 


1 < \k\ < n- 1 


1 


|fc| = n 


1 

-1 



4 



The quantity u>k is defined by 



COS OJk 



1 



2\k\ 



(9) 



Notice that the eigenvalues depend only on n and on the 
Hamming weight of k, defined as |fc| = Xw=o ^j- 

One may define the velocity of a classical walker as the 
derivative of the Hamming distance as function of time. 
We define accordingly the group velocity of the quantum 
walker as 



dto 
d\k\ 



(10) 



where u) is the angular frequency. By examining the 
eigenvalue table, we see that the group velocity is not 
zero only if < |fc| < n, and it is given by 



±1 



y/\k\(n-\k\) 



(11) 



Fig. [T] depicts v g as function of the wave number when 
n = 100. The maximum velocity is 1/y/n — 1 when 
\k\ = 1 or \k\ = n — 1 and the minimum is 2/n when 
\k\ = n/2. Fig. [1] also depicts the phase velocity and 
the dispersion relation. For k < 85 the phase velocity is 
greater than the group velocity. For small k the disper- 
sion relation has negative concavity and for 15 < k < 85 
it is close to a straight line. Those facts indicate that the 
maximal group velocity which is achieved at k = 1 is the 
signal velocity. Then the time for the walker to go to the 
opposite corner of the hypercube is n/v™ ax ~ n^Jn when 
n is large. It is interesting to compare this time with the 
current definitions of hitting time. 



■ group velocity 
- dispersion relation 



phase velocity 




Figure 1: Group velocity, phase velocity and uik of the quan- 
tum walk on the hypercube with n — 100 as function of |fc|. 
The axes are in log scale. 

Taking 4= Y^=o \ e j) ® 1 0, • • • ,0) as the initial condi- 
tion, the one-shot hitting time from vertex (0, • • • ,0) to 



(1, • • • , 1) is either [f n\ or [f n] for p = 1-0(!2S_J2) with 
the condition that the hitting time and n have the same 
parity Q. Note that, as discussed in Q, one cannot in- 
crease the threshold probability beyond p = 1 — 0( og " ) 
without getting infinite one-shot hitting times, because 
the threshold probability is very close to the maximum 
value of the probability distribution at the final vertex. 

Now when n is large, with high probability the walker 
hits the opposite corner at time 0(n). This is faster than 
the 0{n^Jn) scaling of the group-velocity hitting time. 
We will discuss the interesting origin of the difference in 
the scaling behavior in the last section. 

Also, using the same initial condition given above, one 
obtains that the walk has Ct( . 1 — )) concurrent 

v 2 ' v n log n'' 

hitting time [7[. Note that in this case the probability 
of finding the walker at the final vertex is close to zero 
for large n. This result is not a contradiction with the 
group-velocity hitting time for the unitary walk because 
the evolution in this case is non-unitary, i.e., the walk 
is a different physical process, due to the repeated mea- 
surements demanded by the definition of the concurrent 
hitting time. 

In comparison, also the value of the average hitting 
time obtained in Ref. T(j is smaller than the group- 
velocity hitting time. This again is not a contradiction, 
because the approach of [10( does also not describe the 
same physics as we do here, due to the non-unitary evo- 
lution caused by the repeated measurements assumed in 
the definition of the average hitting time [l(| • 



V. QUANTUM WALKS ON A 1-D LATTICE 

In this example, the group-velocity based hitting time 
and the one-shot hitting time are essentially in agreement 
for a suitable choice of the threshold probability. 

An one-dimensional lattice is the Cayley graph of the 
additive group of integers Z with S = {1,-1} as the 
generating set. Since Z is infinite, the theory of Sec. |TT] 
does not apply straightforwardly in this case. We make 
the necessary modifications in this Section. 

The shift operator reduces to 



S\j)\n) = \j)\n + j) 



(12) 



Without loss of generality, we can use the Hadamard ma- 
trix as coin operator (Tl| . which is given by 



C = 



1/1 1 



y/2 \ 1 -1 J ' 

The generic state of the walk is given by 

M*)>= E f>i,»(*)li>l">, 

J={1,— 1} n=-°° 

and the probability distribution by 

p n (t)= J2 i^»wi 2 - 

i={i -i} 



(13) 



(14) 



(15) 
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The transformed amplitudes are 

oo 



(16) 



where k 6 [— 7r, it]. The reduced evolution operator 
which acts on |^ fe ) = £V ^j,feU), is given by 



(17) 




In the eigenbasis, ?7fc is given by 




(18) 



where = e JWfc and 



the angle in [— 7r/2,7r/2] such that sin(u;fc) = sin(fc)/v / 2- 
The Hamiltonian associated with that evolution oper- 
ator is 




Hh, — 



We can now calculate the group velocity v g — 4^ 
w is the angular frequency, obtaining 



± ±cos(fc) 



-y/l + COS 2 (fc) 



(19) 



where 



(20) 



The phase velocity are characterized by two values, which 
are 



1 



. sin k . 



v P h = 1 arcsin(-^-) 



(21) 



and v P h = - v ph - n / k - 

Fig. [5] depicts w+ as function of the wave number. The 

maximum velocity is 1/ 1 \[2 when k = and the minimum 
velocity is the opposite value when k = ±7r. The phase 
velocity (Vp h ) is equal to the group velocity for k = 
and is greater than the group velocity when k > 0. The 
dispersion relation has negative concavity for < k < ir. 
Those facts indicate that the group velocity is in fact 
the signal velocity. For the second values of phase and 
group velocities, we have v~ h < i>~ when < k < n. 
This anomalous case involves velocities that are smaller 
or equal to the maximum group velocity at k = 0. 

It is interesting to relate the group velocity with the 
probability distribution. Fig. [3] depicts the probability 
distribution of the Hadamard walk at t — 100. Note 
that the distribution is clearly non-zero in the region 
—v™ ax t < n < v™ ax t. One can verify that the prob- 
ability distribution is not exactly zero for |n| > v™ ax t 
although very small. 

It is trivial to calculate the one-shot hitting time for the 
lattice case. One can read it directly from the probability 



■ group velocity phase velocity 

- dispersion relation 




Figure 2: Group velocity, phase velocity and uou of the quan- 
tum walk on the lattice as function of k. 




-100 -80 -60 -40 -20 20 40 60 80 100 
II 

Figure 3: Probability distribution of the Hadamard walk at 
t = 100 with initial condition |^(0)) = ( |0> ^ |:L> ) \n = 0). 



distribution. For example, in Fig. [3] we see that the plot 
has a sharp peak at around n max — t/y/2. If we take p 
as the value Pn malL (t) obtained from Eq. (|15p . which is 
the most natural one to take, the one-shot hitting time 
is \/2n. In this case, the one-shot and the group- velocity 
hitting time yield the same value approximately. 

The calculation of the average hitting time is somewhat 
tricky, because it is defined in terms of an only slowly 
converging series. Our numerical results indicate that the 
average hitting time is again smaller than group- velocity 
hitting time. 

Note that parameter n for the line has a different mean- 
ing when compared to the parameter n for the hypercube. 
In the line, n is a linear distance to the origin while in 
the hypercube n is a dimension. 
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VI. DISCUSSION AND CONCLUSIONS 

Quantum walkers often behave similarly to wave pack- 
ets in media. It is the case, in particular, when the 
walker's lattice possesses an Abelian symmetry group 
which allows the use of a normal mode decomposition. 

In this situation, the walker can be described as a wave 
packet which over time propagates and disperses. When 
the walk is unitary the walker's wave packet effectively 
travels in a medium which is neither absorptive nor am- 
plifying. For such passive media, the group velocity is 
known to be a good measure of the speed with which the 
wave packet can propagate information. 

This motivated us to use the group velocity as the ba- 
sis for a new definition of hitting time. The new hitting 
time is defined as the distance divided by the maximal 
group velocity for any wave number, i.e., also for any 
wave packet. Therefore, among all possible initial condi- 
tions for the quantum walker's wave packet, the group- 
velocity based definition of hitting time yields the opti- 
mum. The group velocity based hitting time also does 
not depend on a choice of threshold probability. Instead, 
this hitting time depends only on the coin operator and 
on the symmetry group which defines the Cayley graph. 

Within this approach, we calculated the hitting times 
for discrete-time quantum walks on Cayley graphs of gen- 
eral Abelian groups. For the special cases of the hy- 
percubes and the one-dimensional lattice, we compared 
the group velocity based hitting times with hitting times 
obtained with respect to previous definitions of the hit- 
ting time. While we found general agreement in the case 
of the quantum walk on the line, we found for the hy- 
percubes that the group velocity based hitting times are 
generally scaling slower with n (the dimension) than the 
hitting times with respect to previous definitions. 

To explain the apparent discrepancy, let us first con- 
sider the fact that the group velocity based hitting time 
is larger than the concurrent and average hitting times. 
That there is a discrepancy is not surprising, because the 
described physical processes are different. In the case 
of the group velocity based hitting time calculation, the 
quantum walk is unitary while in the other cases the 
quantum walk is non-unitary due to the performance of 
measurements. 

More significant and interesting is the fact that, for 
the same unitary quantum walk on the hypercube, the 
group velocity based hitting time is larger than the one- 
shot hitting time. 

The group velocity based hitting time is determined 
by how fast the fastest wave packet can travel. In com- 
parison, the one-shot hitting time is based on the idea 
that the arrival of the walker can be recognized by the 
arrival of a certain threshold probability. In the case 
studied in the literature, where the walker is asked to 
reach the diagonally opposite vertex in the hypercube, 
the threshold probability was optimized and it is in fact 
close to one. It would appear, therefore, that the arrival 
of such a large threshold probability indicates the arrival 



of a wave packet. How, therefore, can the group velocity 
based hitting time scale slower than the one-shot hitting 
time? 

To this end, let us consider the larger picture. In prin- 
ciple, in eventual practical applications in quantum com- 
puting, it is important how fast the walker can arrive at 
any vertex - not only the vertex opposite to the starting 
vertex. What, however, is the one-shot hitting time with 
respect to the walker's arrival at a vertex other than the 
one diagonally opposite? The calculation of the one-shot 
hitting time for arrival at the diagonally opposite ver- 
tex gives a partial answer. It was shown there that the 
threshold for arrival can be chosen very high, in fact con- 
verging to 1. The walker is exceedingly likely to arrive 
there. This also shows that the threshold for the walker's 
arrival at other vertices must be chosen small in order to 
obtain a finite value for the one-shot hitting time. In the 
larger picture, where we ask how fast a quantum walker 
can visit any vertex in the graph, this indicates that the 
one-shot hitting time is a difficult measure to use. This 
is because suitable threshold probabilities and therefore 
the one-shot hitting times can heavily depend on the end 
vertex considered. The reason for that is apparently the 
possibility of strong destructive or constructive interfer- 
ences, that lead, in particular, to a very significant en- 
hancement of the arrival probability at the diagonally 
opposite vertex. 

The group velocity based hitting time, on the other 
hand, does not require the consideration of threshold 
probabilities. Nor does it seem to depend on whether or 
not the initial and final positions of a wave packet that 
represent the walker are in a highly symmetric relation- 
ship such as being diagonally opposed to another. In the 
larger picture, where we ask how fast a quantum walker 
can visit any arbitrary vertex of the graph, the group ve- 
locity based hitting time should therefore provide a more 
reliable measure of that speed. 

Nevertheless, for completeness, let us address the re- 
maining question regarding the special case of the walker 
arriving at the diagonally opposite vertex. How can it be 
that, as the one-shot hitting time calculation has shown, 
the walker can quite reliably arrive at the diagonally op- 
posite vertex faster than the group velocity would indi- 
cate is possible? 

To see this, let us recall that wave packets tend to pos- 
sess leading small amplitude waves that travel faster than 
the group velocity. These are the Brillouin and Sommcr- 
feld precursors, which can also be viewed as evanescent 
waves. In general, in passive media, in normal circum- 
stances, precursors stay small during the propagation. In 
active media, precursors may be amplified, thereby lead- 
ing to an apparent speed up of the wave packet. We con- 
jecture that a similar situation prevails here, even though 
the medium is passive (since the evolution is unitary). 
Namely, it could be that the precursors of the quantum 
walker's wave packet, on its way from the initial vertex 
to the diagonally opposite vertex, constructively interfere 
so as to lead to an effectively amplified precursor arriv- 
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ing at the diagonally opposite vertex, before the arrival 
of the main wave packet. This, and also the relationship 
between the group velocity and the notion of mixing time 
should be interesting to explore further. We remark that 
an earlier conjecture by one of us (AK) has been con- 
firmed, see [12J, that the precursor phenomenon occurs 
for the quantum walk on the line. 

It should also be interesting to determine the group 
velocities for more general quantum walks. The method 
we have described here is applicable to Cayley graphs of 
Abelian groups. A suitable generalization may be appli- 
cable to non-Abelian Cayley graphs, or to even to more 
general graphs that allow quantum walks, as long as these 



possess some form of normal mode decomposition. 
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